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ABSTRACT 

We present the Submillimeter Array (SMA) 157-pointing mosaic in 0.86 mm dust continuum emis- 
sion with 5 // .lx4 // .2 angular resolution, and the NRAO Green Bank 100m Telescope (GBT) observa- 
tions of the CS/C 34 S/ 13 CS 1-0 and SiO 1-0 emission with <20"xl8" angular resolution. The dust 
continuum image marginally resolves at least several tens of 10- 10 2 M dense clumps in the 5' field 
including the circumnuclear disk (CND) and the exterior gas streamers. There is very good agreement 
between the high resolution dust continuum map of the CND and all previous molecular line obser- 
vations. As the dust emission is the most reliable optically thin tracer of the mass, free from most 
chemical and excitation effects, we demonstrate the reality of the abundant localized structures within 
the CND, and their connection to external gas structures. From the spectral line data, the velocity 
dispersions of the dense clumps and their parent molecular clouds are ~10-20 times higher than their 
virial velocity dispersions. This supports the idea that the CND and its immediate environment may 
not be stationary or stable structures. Some of the dense gas clumps are associated with 22 GHz 
water masers and 36.2 GHz and 44.1 GHz CH3OH masers. However, we do not find clumps which 
are bound by the gravity of the enclosed molecular gas. Hence, the CH3OH or H2O maser emission 
may be due to strong (proto) stellar feedback, which may be dispersing some of the gas clumps. 
Subject headings: Galaxy: center — Galaxy: structure — Galaxy: kinematics and dynamics — ISM: 
clouds 



1. INTRODUCTION 

The Galactic center (see Morris & Serabyn 1996 and 
Mezger et al. 1996 for reviews) is a fertile ground for 
studying the interplay between the supermassive black 
hole (SMBH; see Genzel et al. 1997 and Ghez et al. 
2005) and the surrounding interstellar medium (ISM). 
The overall dynamics, the dumpiness, and the local kine- 
matics may provide clues as to how the SMBH is fed by 
the ISM, how the circumnuclear gas streams evolve, and 
how the molecular cores and the young massive-stellar 
objects (Krabbe et al. 1991; Launhardt et al. 2002; Pfuhl 
et al. 2011, Lu et al. 2013, Do et al. 2013, etc) form 
and migrate in the Galactic center. Previous molecular 
line observations indicated that the warm and turbulent 
gas clouds or streamers in the central ~20 pc in Galactic 
center (Ho et al. 1991; Okumura 1991; Coil & Ho 1999, 
2000; McGary et al. 2001; Herrnstein & Ho 2002, 2005; 
Oka et al. 2011) are well connected with the 2-4 pc cir- 
cumnuclear disk (CND; Glisten et al. 1987a; Jackson et 
al. 1993; Marshall et al. 1995; Christopher et al. 2005; 
Montero-Castaho et al. 2009; Liu et al. 2012; Martin et 
al. 2012) surrounding the central black hole. These gas 
structures are the most extreme environments for high- 
mass star formation in the Milky Way (Morris 1993 and 
references therein). 

The thermal dust continuum emission is the least bi- 
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ased by the chemical abundances and the excitation con- 
ditions. Hence, it is the most reliable tracer for demon- 
strating the reality and significance of the individual gas 
clouds or streamers, as well as their internal structures. 
The IRAM-30m and the JCMT single dish telescope ob- 
servations of the millimeter and the submillimeter con- 
tinuum emission have been presented by Zylka & Mezger 
(1988), Mezger et al. (1989), Dent et al. (1993), Lis 
& Carlstrom (1994), Pierce-Price et al. (2000), and 
Garcia-Marm et al. (2011). In this work, we report the 
first wide-field inter ferometric mosaic observations (157 
pointings, ~5'x5' field of view) of the 0.86 mm dust con- 
tinuum emission using the Submillimeter Array (SMA0; 
Ho et al. 2004). The improved angular resolution of ~5" 
permits the detection of 0.1-0.2 parsec scale gas clumps, 
which can be the candidates of high-mass star-forming 
cores. In addition, we compare the dust continuum im- 
age with the HCO + 4-3 line image simultaneously ob- 
tained with the SMA. We also compare these maps with 
the National Radio Astronomy Observatory (NRAO0) 
Robert C. Byrd Green Bank Telescope (GBT) observa- 

1 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 

2 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



tions of the CS/C 34 S/ 13 CS 1-0 lines and the SiO 1-0 lines. 
This allows us to examine the gravitational stabilities of 
the dense structures and the possible existence of shock 
fronts (e.g. Kauffmann et al. 2013). 

The new observations are described in Section [2j The 
results are presented in Section [U We compare our obser- 
vations with the previous observations of OH, CH3OH, 
and H2O masers in Section |U A discussion of the non 
uniformity of the CND based on the 0.86 mm dust con- 
tinuum image is given in Section 14.41 A brief summary 
of our results is given in Section 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Submillimeter Continuum Emission 

2.1.1. SMA Observations 

We made mosaic observations toward the Galactic cen- 
ter using the SMA, in its compact and subcompact array 
configurations. These observations covered the frequency 
range of 354.1-358.1 GHz in the upper sideband, and 
342.1-346.1 GHz in the lower sideband. Details of the 
subcompact array observations, and the pointing centers 
of the observations, can be found in Liu et al. (2012). 
The compact array observations were made in two ob- 
serving runs, on 2012 May 07 and 2012 May 20, with 6 
and 7 available antennas, respectively. The system tem- 
peratures T sys during these two runs were ~ 180-400 K. 
We observed two phase calibrators, 1733-130 and 1924- 
292, every ~15 minutes during all observations. The 
amplitude and passband calibrators were Titan and 3C 
279 on May 07, Neptune and 3C 279 on May 20. The 
target loops iterated over 157 pointing centers, with 5 
5-second integrations at each pointing center. Each of 
the 157 pointings was visited more than three times in 
each observing run (i.e. on source time >3 hours in 
total in each run). The minimum and maximum pro- 
jected baselines in our SMA observations are ~7.0 kX and 
-82 kX. All SMA data were calibrated using the MIR 
IDL software package (Qi 2003). We used the MIRIAD 
(Sault 1995) task UVAVER to average all line- free chan- 
nel data and reconstruct the 0.86 mm continuum band 
data. The short spacing data were complemented by 
combining the archival JCMT SCUBA image (Appendix 
A, B). The zeroth order free- free continuum emission 
model was constructed based on the archival VLA 7 mm 
observation data, and was subtracted from the combined 
SMA+JCMT 0.86 mm continuum image (Appendix C). 
The simultaneously observed HCO + 4-3 line was regrided 
to 2.1 krns" 1 velocity channels for an adequate sensitiv- 
ity, and will be presented for the purpose of discussing 
the virial condition. The detailed studies of the submil- 
limeter line data are deferred. 

2.2. GBT Observations 

We observed the CS 1-0 (48.99095 GHz) and the C 34 S 
1-0 (48.20692 GHz) transitions using the NRAO GBT 
on 2011 November 04 and 07. We observed the 13 CS 1-0 
(46.24754 GHz) and the SiO 1-0 (43.42376 GHz) transi- 
tions on 2011 November 07 and 09. The field of view of 
the SiO and 13 CS observations is slightly offset toward 
the west to better recover the western gas streamers. 
The angular resolution of the GBT is 763.8"/^, where 
v is the observing frequency in GHz. The bright point 
source 1733-130 was observed in the beginning of each 



session for antenna surface (i.e. by Out of Focus Holog- 
raphy) and pointing calibrations. A line- free reference 

position RA: 17 h 43 m 43 s .344, DecL: -29°59 / 32 // .27 was 
integrated for 30 seconds before and after the target ob- 
servations in each block for off-source calibration data. 
We used the GBTIDL software package (Marganian et 
al. 2006) to calibrate the GBT data. We note that the 
CS isotopologue lines are close to the band edge of the 
GBT, which are subjected to the higher system temper- 
ature. We used the AIPS software package to perform 
imaging. We smoothed the final CS and C 34 S image to 
an optimized ma jX0 m i n = 20 // xl8 // , and BPA. = 0° to 
suppress the striping defects due to the sampling rates, 
and removed some obvious stripes by fitting zeroth order 
polynomial. The achieved RMS noise levels in each 24 
kHz (~0.15 kms -1 at 48.99 GHz) spectral channel are 
~0.5 K for the CS and C 34 S observations, and are ~0.14 
K for the 13 CS and SiO observations. 

3. RESULTS 

We present the observing results in this section. Our 
nomenclature follows Christopher et al. (2005), Amo- 
Baladron et al. (2011), and Liu et al. (2012). 

We also compare our results with the VLA observa- 
tions of the 1612 MHz OH masers (Pihlstrom et al. 2008; 
see also Sjouwerman 1998), the VLA observations of the 
1720 MHz OH masers (Sjouwerman & Pihlstrom 2008), 
the GBT and VLA observations of the 44.1 GHz Class I 
CH3OH masers (Yusef-Zadeh 2008), the JVLA observa- 
tions of the 44.1 GHz Class I CH3OH masers (Pihlstrom 
et al. 2011), the JVLA observations of the 36.2 GHz 
Class I CH3OH masers (Sjouwerman 2010), and the GBT 
observations of the 22 GHz H2O masers (Yusef-Zadeh 
2008). 

3.1. Continuum Emission 

Figure Q] shows the high angular resolution 
SMA+JCMT 0.86 mm continuum image. A model 
of the free-free emission has been subtracted (see 
Appendix IC]) . The 0.86 mm continuum image recovers 
the detailed gas structures in the CND including the 
Northeast Lobe, the Northeast Extension, the South- 
west Lobe, the Southern Extension, and the W-2,3,4 
Streamers which connect to the CND from the west. It 
addition, it presents the detailed structures embedded 
in the Northern Ridge, the Southern Ridge, the 50 
kms -1 cloud, the Molecular Ridge, the 20 kms -1 
cloud. A southern dust ridge which appears to connect 
JCMT-1 to JCMT-2 (see Figure CD ) in the archival 
JCMT SCUBA 0.44 mm (678 GHz) continuum image 
published by Pierce-Price et al. (2000; ProjectID: 
M98AU64), is also resolved in our SMA+JCMT image 
(see Figure [2]). The spatial location of this southern 
dust ridge coincides with the Southern Arc reported by 
the previous observations of the CS 1-0 line emission 
(Liu et al. 2012). 

Excluding the central point source, the 0.86 mm flux 
of the CND is approximately 68 Jy in a distribution that 
closely follows that of the inner 5 pc of the CND. At least 
few tens of dense molecular gas clumps are marginally 
resolved. Higher angular resolution observations may 
resolve more blended dense clumps in this field. The 
precise number of clumps is not important for this dis- 
cussion. The main conclusion is the impression of a very 
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Fig. 1. — The SMA+JCMT 0.86 mm continuum image with a free-free model subtracted, (color and contour). The angular resolution 



of this image is n 



5 // .lx4 // .2. Contour spacings are 3a starting at 3a (a=24 mJybeam 1 ). Yellow crosses are the 1720 MHz 
Pihlstrom (2008). Yellow triangles are the compact, either thermal or low-gain masing 1612 MHz 



OH masers taken from Sjouwerman 

OH line sources discussed in Pihlstrom et al. (2008) (see also Sjouwerman 1998). Pink diamonds are the 36.2 GHz Class I CH3OH masers 
reported in Sjouwerman (2010). Yellow diamonds are the 44.1 GHz Class I CH3OH masers reported in Pihlstrom et al. (2011). White 
Crosses and diamonds are the 22 GHz water masers and the 44.1 GHz CH3OH masers reported in Yusef-Zadeh (2008). Note the symbol 
size is not representative to the spatial uncertainty. For our convenience in discussion, we mark the peaks above the 45a significance level 
besides the Northeast Lobe and Southwest Lobe as JCMT-1,2,3 (see Appendix A). The regions associated with CH3OH masers are labeled 
by A-H. 



clumpy structure. The most significant dense clumps 
are found in the protrusion connecting to the Southwest 
Lobe. Clusters of very dense clumps are also resolved 
towards the peaks of the JCMT SCUBA 0.86 mm image 
(e.g. JCMT-1,2,3 and the peaks in the Northern Ridge; 
see Appendix A). 
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Fig. 2. The JCMT SCUBA 0.44 mm (678 GHz) continuum 
image (color; # ma j x m ^ n =8 n x8 n ), overlaid with the high angu- 
lar resolution SMA+JCMT 0.86 mm continuum image (contour; 
^maj x ^min— 5 // .lx4 // .2). Contour spacings are 3cr starting at 
3cr (cr=72 mJybeam- 1 ). The public processed JCMT SCUBA 
0.44 mm was retrieved from the online data archive (Project ID: 
M98AU64), and was published by Pierce-Price et al. (2000). The 
beam of the JCMT SCUBA 0.44 mm observation is shown in the 
lower left. 



In the present work, only the dense clumps located near 
the CH3OH masers, and additionally the dense clumps 
located in JCMT-1 and JCMT-2 are discussed (Figure 
[3|). We do not attempt to systematically search for and 
analyze all dense clumps in the entire field because of 
the limited image quality, and insufficient angular res- 
olution for the CND (see Montero-Castaho et al. 2009 
and Martin et al. 2012 for higher angular resolution spec- 
tral line images). To estimate the 0.86 mm fluxes of the 
selected clumps, we first fit ~1 pc scale 2D Gaussian 
components to suppress the contribution of the ambi- 
ent gas, and then fit 2D Gaussian components to the 
localized peaks. The 0.86 mm fluxes of the individual 
Gaussian components are summarized in Table [TJ which 
provides a quantitative measure of how massive the sub- 
parsec scale gas over-densities may be. The uncertain- 
ties in background subtraction can give a few tens of 
percent of errors in fluxes. In addition, the overesti- 
mates of the background emission can lead to underes- 



timates of the clump sizescale. The main purpose of 
this analysis is to estimate the virial velocity dispersions 
of these dense clumps (Section I4.3J) . In this sense, the 
effects of overestimates of the background and underes- 
timates of the clumps sizescale are competing. Gaussian 
components with minor axis FWHM smaller than one 
standard deviation of the SMA+JCMT Gaussian beam 
width, 5". 2/2. 355 = 2.2" are not considered because we 
cannot estimate the physical deconvolved size scales. The 
Gaussian components which are located near the edge of 
the SMA+JCMT 0.86 mm continuum image are also not 
considered. 

3.2. Molecular Emission 

The SMA image of HCO+ 4-3 is shown in Figure gj 
over-plotted with the fitted 0.86 mm dense clumps (Sec- 
tionOTTJ Figure[3]). The HCO+ 4-3 line is a good tracer of 
the ring-like CND, and also the gas streamers connecting 
to the CND. The eastern part of the CND is relatively 
narrow and smooth, compared to the western part. The 
CND may be undergoing dynamical evolution, or is com- 
posed of distinct streams of molecular gas (Jackson et al. 
1993; Wright et al. 2001; see Section WM- 

Towards the more extended 20 kms -1 cloud, the 
Southern Arc, the Molecular Ridge and the 50 kms -1 
cloud, the SMA observations are detecting strong miss- 
ing flux (see Liu et al. 2012). In addition, the lower gas 
temperature in the extended cloud than in the CND also 
causes the weak or non-detection of the warm gas tracer 
HCO+ 4-3 (E up = 43 K). We therefore can only robustly 
image the HCO + 4-3 line towards a few fitted dense 
clumps. The detected HCO + 4-3 spectra are presented in 
Figure [5] The velocity dispersions derived from the sin- 
gle component Gaussian fittings are summarized in Table 
[H From Figure [3] and [5] we see that the brightest clump 
over the selected samples (Table [1]), the JCMT-ln, has 
a centroid velocity of -95 kms -1 . The JCMT-ln clump 
is likely to be embedded in the very blueshifted South- 
ern Ridge, which shows strong SiO 2-1 emission in the 
earlier Nobeyama Millimetre Array observations (Sato 
& Tsuboi 2008, and also see Amo-Baladron et al. 2011). 
However, the CS 1-0 line spectrum taken at JCMT-ln 
indicates that a vi sr ~7 kms -1 gas component is com- 
pletely missed from the SMA observations of HCO + 4-3 
(ref. spectrum 10 in Liu et al. 2012). The virial velocity 
of JCMT-ln derived from the 0.86 mm dust continuum 
emission (Table [1]) then should be considered as an up- 
per limit when comparing with the HCO + 4-3 velocity 
dispersion. 

The more extended gas streamers are recovered by the 
GBT observations of the CS/C 34 S/ 13 CS 1-0 lines and 
the SiO 1-0 line (Figure [6]), which trace cooler gas (E up = 
2.1-2.3 K). The SiO 1-0 emission appears to be correlated 
with the CS 1-0 emission. The flux ratio of these lines 
will be discussed in Section 14.11 

4. DISCUSSION 

We discuss the inferred ISM properties based on the 
presented observations. 

4.1. The Spectral Line Ratio 

While the abundance of the CS molecule is only 
mildly enhanced in UV and shocked environments (Amo- 
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Fig. 3. — Blowups of the free-free model subtracted SMA+JCMT 0.86 mm continuum image (color and contour) in parsec scale area 
around regions A, B, D, E, F, G, JCMT-1, and JCMT-2. Contour spacings are 3cr starting at 3cr (cr=24 mJybeam -1 ). Red ellipses show 
the fitted 2D Gaussian clumps. The major and minor axes of the red ellipses are the FWHM of the fitted 2D Gaussian (not deconvolved). 



Baladron et al. 2011 and references therein), the abun- 
dance of the SiO molecule in shocked environments can 
be enhanced by up to a factor of 10 6 with respect to 
the value in quiescent gas (Mart in- Pint ado et al. 1992). 
Since SiO and the CS isotopologues have similar dipole 
moments and energy level distributions, the derived flux 
ratio between these molecules can trace shocks without 
being sensitive to the assumed physical conditions. 

Figure shows the velocity integrated CS/C 34 S 1-0 
line ratio. Above ~40 GHz, difficulties in calibrating the 
several arcsecond GBT pointing offsets can potentially 
lead to the ~30% uncertainty in absolute flux level. This 
uncertainty does not affect the derived CS/C 34 S 1-0 ratio 
since these two lines are simultaneously observed. We 



found that the CS/C 34 S 1-0 ratio is ~6-12 towards the 
50 kms -1 cloud and the 20 kms -1 cloud. The earlier 
NRO 45m Telescope observations of CS 1-0 and C 34 S 
1-0 only robustly detected the C 34 S 1-0 emission at the 
reference point located 3' north and 3' east of the Sgr 
A*, and showed a CS/C 34 S 1-0 ratio of 8 (Tsuboi et 
al. 1999). Considering that the NRO 45m Telescope 
observations were more beam smoothed, our results are 
consistent with the previous observations. Assuming the 
abundance ratio of [X(CS)/X(C 34 S)] = 22.6 (Frerking 
1980; Tsuboi et al. 1999), the optical depth r of CS 1-0 
can be estimated based on the following relation 
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Fig. 4. — The velocity integrated HCO+ 4-3 image (<9 ma j x 
$min— 5 // .8x4 // .l). Red ellipses show the fitted 0.86 mm clumps 
(see Figure [3}. Blue star labels the location of Sgr A*. 



10 
8 
6 

4 
2 

0.0 
10 
8 
6 
4 
2 

0.0 
10 
8 
6 
4 
2 

0.0 
10 
8 
6 
4 
2 

0.0 
T 10 

S 6 

\ 4 

3 2 
- 0.0 



Region A 



ry^uu^-^w^ n^^n^^^ 




Region Bn 



- — q JLJ1 ^ J ^rv H ^^^n^r JWx nj^ri 1 ^ - LJ : jij^^ 




rJl^j . „ J ^j^u LnJ~^Il 



Region B 



^ hJ~ ~iQ^ n -H~rJ~Hx^^ 



-^^V 




■ "Tr-^piJ 1 - ^-njxrJ^T-T r 



Region Bs 



^^^^U~^^^^ ^ T_Jlj ° ^ 



_ JCMT-1n 



r ~ J lj ^ Lr~ J -njTjxr^ 



-150 -100 -50 50 

Velocity (km/s) 



1 00 1 50 



Fig. 5. — The HCO + 4-3 line spectra toward the selected re- 
gions (black; see Table [T}. These spectra are regridded to 2.1 
kms -1 velocity channels. The rms noise level is 0.6 Jybeam -1 (1 
Jybeam _1 ~0.4 K). Blue lines present the fitted Gaussian compo- 
nents. 



where Scs and S C 34 s are the fluxes of the CS and the 
C 34 S lines. The line ratio 5cs/# c 34 s= 14 - 2 when r=l. 
From Figure [3 the majority of the optically thick CS 
1-0 line (r=1.5-4, i.e. CS/C 34 S 1-0 ratio -6-12) is seen 
towards the 50 kms -1 cloud and the 20 kms -1 cloud. 
The C 34 S 1-0 line is optically thin over the observed field. 

The SiO 1-0 line is —1.5-2 times fainter than the CS 
1-0 line (Figure [6j). Assuming the optically thin limit, 
the velocity integrated SiO/C 34 S 1-0 line ratio (Figure 
[8]) traces the abundance ratio. We present the velocity 
channel maps of the SiO/C 34 S 1-0 line ratio in Figure [9] 
to trace the velocity structures of the shocked gas. While 
the SiO/C 34 S 1-0 ratio is already high in the entire map 
area, it appears to be enhanced in the 20 kms -1 cloud 
and the Southern Arc. Based on the correlation between 
NH3 and 1.2 mm dust emission and the locations of OH 
masers, Wright et al. (2001) also suggested that these 
regions could be shock heated. Our observations do not 
have sufficient angular resolution to resolve the shock 
fronts. Whether the shocks are created by the super- 
nova shell to the south of the Sgr A East, or are created 
due to the non-circular orbits of the clouds, are not yet 
distinguished. 

The earlier IRAM-30 m Telescope observations addi- 
tionally showed the enhanced SiO emission toward the 
50 kms -1 cloud and around the Southern Ridge and the 
W-4 Streamer. Our GBT observations of SiO 1-0 only 
covered a small part of the 50 kms -1 cloud, and did 
not significantly detect the Southern Ridge. The W-4 
Streamer is only marginally detected in our C 34 S 1-0 ob- 
servations. 

4.2. The Molecular Gas Environment Around the 
Masers 

The 1720 MHz OH masers and the 1612 MHz OH 
masers trace the nn 2 — 10 5 cm -3 and nn 2 >10 7 cm -3 
postshock gas in supernova remnants (assuming T— 75 
K; Pihlstrom et al. 2008; Pavlakis & Kylafis 1996). The 
22 GHz H2O masers are collisionally excited at higher 
densities (nn 2 ~10 7-9 cm -3 ; Elitzur et al. 1992), and 
are often found in star- forming regions. The collisionally 
excited Class I CH3OH masers are regarded as the unam- 
biguous signposts of ongoing high-mass star formation in 
typical molecular clouds (Menten 1992; Yusef-Zadeh et 
al. 2008). However, the Galactic center (e.g. inner 2-20 
pc) gas streamers are warmer and 10 times more tur- 
bulent than the other star- forming molecular clouds. In 
such extreme environments, the Class I CH3OH masers 
may also be excited in cloud interactions although not 
yet observationally confirmed (Sjouwerman et al. 2010). 

From Figure [T3l we can see that on —10 pc scales, 
while the distribution of the OH masers follows a shell 
like geometry, the CH3OH and H2O masers are more 
scattered over the field. In the higher angular resolu- 
tion dust continuum image (Figured]), we see that the 
CH3OH and H2O masers are preferentially detected near 
the <0.5 pc scale localized over-densities. The only ex- 
ceptions are the CH3OH masers associated with the re- 
gion H which is located near the central bright free-free 
continuum sources. The flux at region H may be biased 
by the uncertainties in interferometric imaging and the 
free- free model subtraction (Appendix IC]). In addition, 
the molecular gas in region H may be photo-ionized by 
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Fig. 6. — Molecular gas in the central 20 pc region in the Galactic center. (A) The velocity integrated CS 1-0 emission. (a ~0.7 Kkms 1 ). 
The symbols are described in Figure [T] Labels of the large scale gas streamers are consistent with those in Liu et al. (2012). (B) The 
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13 CS peaks at R.A.=17 45 m 34-35 s may be caused by ambiguities in baseline subtraction. The rms noise level of the velocity integrated 
images are estimated based on the integration of the signal over a 13 kms -1 velocity range, which is the median of the fitted velocity 
dispersion (see Table ^). Green contours in panels (B) and (D) show the free-free model subtracted SMA+JCMT 0.86 mm continuum 
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Fig. 7. — Velocity integrated CS/C 34 S 1-0 ratio (color and con- 
tour). Contour levels are [6, 12, 18, 24]. Blue star labels the 
location of Sgr A*. 

the central engine Sgr A*. 

We summarize the observed intensity of the dust con- 
tinuum emission and the velocity integrated SiO/C 34 S 
1-0 ratio in Figure [p Except for the 44.1 GHz CH 3 OH 
masers associated with region H, the other CH 3 OH and 
H2O masers are detected at intensities >0.15 Jybeam -1 . 
The OH masers do not show a clear correlation with the 
flux of the dust continuum emission. The clusters of 
OH maser may only sparsely sample the geometrically 
thin expanding shell (Coil & Ho 1999) in projection of 
the bulk of the dense gas. Yusef-Zadeh et al.(2008) and 
Pihlstrom et al.(2011) only observed H 2 and CH 3 OH 
masers towards selected high density regions, which po- 
tentially bias the corresponding maser data in Figure [TOl 
toward higher 'averaged' 0.86 mm intensity. Neverthe- 
less, the size scales of the dense clumps in Figure [1] are 
generally smaller than what can be discerned by previ- 
ous observations of molecular gas. We therefore do not 
think the spatial correlation between the clumps and the 
maser spots is purely an artifact due to the observational 
selection. As an example, the good correlation between 
the 0.86 mm emission clumps in region G (Figure [TJ i.e. 
pointing E in Pihlstrom et al. 2011) and the JVLA de- 
tections of 44.1 GHz CH3OH maser spots does not seem 
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Fig. 8. — Velocity integrated SiO/C 34 S 1-0 ratio (color con- 
tours). The yellow, orange, and red contour levels are 2, 4, 
and 6, respectively. Gray contours show the free- free model sub- 
tracted SMA+JCMT 0.86 mm continuum image. The velocity 
integrated SiO 1-0 emission is shown in grayscale. The 0.86 mm 
continuum image contour spacings are 7.5cr starting at 7.5cr (cr=24 
m Jybeam - 1 ). The GBT images are smoothed to the angular res- 
olution of # ma j x # m i n = 20" x 18" before taking the ratio. Blue star 
labels the location of Sgr A*. 



to be a coincidence, but should be interpreted with the 
local physical conditions. 

For both the maser sources and the ambient gas, our 
observations do not resolve a clear correlation between 
the [X(SiO)/X(C 34 S)] abundance ratio and the 0.86 mm 
intensity (Figure fTO]) . On the large (e.g. >1 pc) scale, the 
shocked gas may be well mixed with the ambient mate- 
rial. The enhanced abundance of SiO in shocks last for a 
few 10 3 years. With the ~50-100 kms -1 relative motion 
of gas clouds around the CND, the smeared local SiO 
shock fronts may have widths of 0.15-0.3 pc (i.e. 3". 8- 
7". 5), which have to be examined with sensitive higher 
angular resolution observations. Here we refer to Martin 
et al. (2012) for higher angular resolution observations 
of SiO emission in the 2' field around the Sgr A*. The 
velocity information is crucial for understanding whether 
and how the shock can induce the formation of sub-parsec 
scale dense gas structures. 
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4.3. The Dense Molecular Clumps 

The brightness temperature of the dust emission at the 
0.86 mm wavelength (Appendix [B]) is much lower than 
the gas temperature (see Herrnstein & Ho 2005). The 
molecular gas mass can therefore be calculated from the 
0.86 mm flux based on the optically thin formula 



M H2 



2X 3 RapD 2 

3hcQ(\)J(\,T d ) 



5(A), 



(2) 



where R is the gas-to-dust mass ratio, a is the mean grain 
radius, p is the mean grain density, D is the distance of 
the target, Q(A)ocA _/3 is the grain emissivity, Td is the 
dust temperature, S(X) is the flux of the dust emission at 
the given wavelength, J(A,T^) = l/[exp(/ic/AfcB^d) — 1] 
(Hildebrand 1983; Lis et al. 1998). The c, h, and k B 
are the light speed, the Planck constant, and the Boltz- 
mann constant, respectively. Following Lis at al. (1998), 
we adopt i?=100, a=0.1 /im, p=3 gem -3 , Q(A=350 
//m)=lxl0 -4 . We adopt D=8.33 kpc based on the mea- 
surements of Gillessen et al. (2009). 

The gas temperature measurements which have a com- 
parable angular resolution with our dust continuum map 
are not yet available. By measuring the rotational tem- 
perature of the NH3 molecule, Herrnstein & Ho (2005) 
suggested that roughly one quarter of the molecular gas 
comprises a hot (~200 K) component, and the remaining 
gas is cool (~25 K). We tentatively assumed an averaged 
gas temperature T^ as =70 K in the dense clump. We also 
assumed the dust temperature Td=T gas (see Chan et al. 



1997 for dust temperature in the CND; see also Minh et 
al.1992, Ao et al. 2012). The ranges of gas mass in the 
dense clump estimated based on the assumption of j3=l-2 
is given in Table [U Based on the aforementioned assump- 
tions, without considering the foreground/background 
subtractions, the detect ed .86 mm flux in the inner 5 pc 
CND (68 Jy, see SectionEl]) corresponds to 0.98-2.3 xlO 4 
Mq of gas mass. Our estimates of the CND mass agree 
reasonably well with those from earlier observations of 
millimeter and submillimeter dust continuum emission 
(e.g. Mezger et al. 1989; see also Christopher et al. 2005 
for relevant debates). 

If the identified gas dense clumps are virialized, the 
expected one-dimensional velocity dispersion Sv 
given by 



viriral :,-, 



Sv 



viriral 



aMG 
5S r 



(3) 



where a is the geometric factor which equals to unity for 
a uniform density profile and 5/3 for an inverse square 
profile (Williams, de Geus & Blitz 1994; Walsh et al. 
2007), M is the gas mass, G is the gravitational constant, 
and 5 r is the effective radius of the dense clump. 

The effective radius S r can be calculated based on the 
FWHM of the fitted 2D Gaussian (i.e. S0 m ^ and S0 min 
in Table [1]), however, is necessary to be corrected for the 
beam FWHM S0 heam because the identified clumps are 
only marginally resolved. For simplicity, we adopt the 
corrected effective circular angular diameter of the dense 
clump to estimate their linear diameter (c.f. Williams, 
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TABLE 1 

Properties of the fitted 22 0.86 mm continuum clumps. 



Components 


R.A. 


Decl. 


se ma i 


se min 


0.86mm Flux 


Gas Mass 


Number Density 


^viriral ^HCO+ ^CS 




(J2000) 


(J2000) 


(arcsec) 


(arcsec) 


(mJy) 


(M ) 


(10 5 cm- 3 ) 


(kms -1 ) (kms -1 ) (kms -1 ) 


Region A 


17 n 45 m 43.74 s 


-29°00 / 03 // 


7.0 


3.2 


90 


13-31 


2.5-5.9 


0.4-0.8 9.7 19 


Region Bn 


17 h 45 m 43.72 s 


-29°00 / 18 // 


5.6 


5.3 


172 


24-60 


1.8-4.4 


0.45-0.96 4 


5 15 


Region B 


17 h 45 m 43.80 s 


-29°00 / 27 // 


6.0 


4.7 


238 


34-83 


2.9-7.2 


0.56-1.1 7 


.0 14 


Region Bs 


17 h 45 m 43.65 s 


-29°00 / 32 // 


5.0 


4.3 


114 


16-40 


4.4-11 


0.46-0.96 9 


4 14 


Region Bwn 


17 h 45 m 42.89 s 


-29°00 / 22 // 


4.1 


3.1 


54 


8-19 






21 


Region Bws 


17 h 45 m 42.81 s 


-29°00 / 28 // 


6.9 


5.7 


317 


45-110 


1.5-3.6 


0.55-1.1 


19 


Region Dn 


17 h 45 m 42.51 s 


-29°01 / 53 // 


9.6 


4.8 


271 


38-94 


0.88-2.1 


0.47-0.96 


14 


Region Ds 


17 h 45 m 42.20 s 


-29°02 / 00 // 


8.2 


4.6 


195 


28-68 


1.0-2.6 


0.43-0.87 


12 


Region Ew 


17 h 45 m 50.28 s 


-29°00 / 06 // 


4.0 


3.0 


121 


17-42 






13 


Region Ee 


17 h 45 m 50.51 s 


-29°00 / 07 // 


10.0 


4.9 


474 


67-160 


1.3-3.2 


0.61-1.2 


13 


Region Fw 


17 h 45 m 49.90 s 


-28°59 / 41 // 


7.1 


3.5 


134 


19-46 


2.6-6.4 


0.45-0.87 


12 


Region Fe 


17 h 45 m 50.16 s 


-28°59 / 44 // 


4.7 


3.5 


83 


12-29 


59-140 


0.64-1.3 


12 


Region Gn 


17 h 45 m 49.41 s 


-28°58 / 54 // 


6.5 


6.5 


364 


51-130 


1.5-3.6 


0.57-1.1 


12 


Region Gs 


17 h 45 m 49.63 s 


-28°59 / 00 // 


4.0 


4.0 


147 


21-51 


320-760 


1.0-2.1 


13 


JCMT-ln 


17 h 45 m 36.56 s 


-29°01 / 58 // 


11.9 


8.0 


1810 


260-630 


1.4-3.4 


0.96-1.9 6 


.8 24 


JCMT-ls 


17 h 45 m 36.10 s 


-29°02 / 07 // 


6.5 


4.9 


527 


74-180 


4.4-11 


0.77-1.6 


29 


JCMT-2a 


17 h 45 m 41.06 s 


-29°02 / 20 // 


5.5 


3.3 


119 


17-41 


16-39 


0.58-1.1 


10 (or 19) 


JCMT-2b 


17 h 45 m 40.22 s 


-29°02 / 24 // 


5.8 


4.2 


169 


24-59 


3.6-8.8 


0.51-1.0 


11 (or 20) 


JCMT-2c 


17 h 45 m 40.83 s 


-29°02 / 26 // 


3.8 


2.5 


48 


6.6-16 






11 (or 18) 


JCMT-2d 


17 h 45 m 40.60 s 


-29°02 / 39 // 


8.8 


4.7 


412 


58-140 


1.8-4.1 


0.61-1.2 


11 (or 17) 


JCMT-2e 


17 h 45 m 41.90 s 


-29°02 / 29 // 


6.8 


5.4 


398 


56-140 


2.2-5.6 


0.63-1.3 


10 (or 23) 


JCMT-2f 


17 h 45 m 39.46 s 


-29°02 / 41 // 


9.4 


5.7 


614 


87-210 


1.4-3.5 


0.68-1.4 


11 (or 26) 



Note.— The 5<9 ma J and 50 min are the FWHM of the fitted 2D Gaussians. The la velocity dispersion 
20" X 18" GBT beam area centered at th e pe ak of the fitted 2D Gaussian. The ranges of gas mass in the 
on the assumption of /3=l-2 (see Section 



Sv is measured in a smoothed 
dense clump are estimated based 



de Geus & Blitz 1994) 



56 ( 



efT 



(56 ma i x S0 min ) - (2£<9 beam /2.355) 2 (4) 



and consider the gas mass enclosed in one FWHM of 
the fitted 2D Gaussian when calculating the virial one- 
dimensional velocity dispersion. For those identified 
dense clumps with 80 e ^ > 0, we compare the results 
of the calculation with the one-dimensional velocity dis- 
persion measured from the CS 1-0 line and the HCO + 
4-3 line (when significantly detected) in Table [TJ The 
mean gas number densities are calculated based on the 
aforementioned estimates of the clump masses and radii, 
and the assumptions of spherical geometry and the mean 
molecular weight of 2.33 (Shull & Beckwith 1982). We 
note that the JCMT-2 region shows double peak pro- 
files in the CS 1-0 spectra. However, the broader line 
component is very faint and therefore is less likely to be 
associated with the dense gas clumps. 

We found that the dense clumps listed in Table [1] have 
~10-10 3 Mq of molecular gas, which can be adequate gas 
reservoirs to form high-mass stars. However, the mean 
gas number density in these dense clumps are generally 
in between 10 5 -10 6 cm -3 , which are marginally unstable 
against the tidal force (Morris 1993; Liu et al. 2012). The 
GBT observations of CS 1-0 suggest that the identified 
dense clumps are all embedded in very turbulent envi- 
ronment, which have one-dimensional velocity dispersion 
of 10-30 kms -1 . The kinetic energy of the gas may be 
dissipated on smaller scales. The higher angular reso- 



lution SMA observations of HCO + 4-3 towards the en- 
closed Region A, Bn, B, Bs, and JCMT-ln detect the 2-3 
times smaller velocity dispersion in the localized dense 
clumps. Nevertheless, the observed HCO + 4-3 velocity 
dispersions are still ~10 times higher than the virial ve- 
locity dispersions of the dense clumps. Our preliminary 
analyses do not yet find any (self-)gravitationally bound 
gas structures on <0.5 pc scale. 

Although our HCO + 4-3 image only traces dense 
clumps in the warmer environments, HCN 1-0 and HCO + 
1-0 spectra from previous 13 // x4 // resolution BIMA ar- 
ray observations also showed broad line profiles in the 
central 5'x5' field (Wright et al. 2001). Observations of 
the lower excitation gas tracers HCN 1-0 and HCO + 1-0 
trace the velocity dispersion of the dense clumps in the 
central 5' 'x5' field (Wright et al. 2001). With 13 // x4 // 
angular resolution, the BIMA spectra are contaminated 
by emission from the diffuse and more turbulent ambient 
gas as traced by CS 1-0, however a deconvolved HCN 1- 
image at ~ 2" resolution (Glisten et al. 1987b shows 
only the integrated image) also shows line widths which 
are much larger than the virial velocity dispersions of the 
dense clumps. Unless these dense clumps are composed 
of dense, bound cores, or the gas kinetic energy can be 
dissipated efficiently, they may be dispersed in one dy- 
namical timescale (i.e. ^ r /Jv~5.3-10 3 -2.6-10 4 years for 
the HCO + 4-3 emission clumps in Table 1). Higher an- 
gular and velocity resolution spectral line observations 
are required to see whether (self-)gravitationally bound 
gas structures exist on a smaller scale. 



11 



Alternatively, these dense clumps embedded in the ex- 
tremely turbulent molecular clouds may be confined by 
the high external pressure. For example, the previous X- 
ray studies (e.g. Koyama et al. 1996, Muno et al. 2004) 
have demonstrated that a high-pressure medium with 
sound speed >1000 kms -1 pervades the Galactic center 
region. If small scale virialized gas structures do exist in- 
side the pressurized dense clumps, then hydrostatic cores 
and stars may still form. Assuming approximate hydro- 
static equilibrium of the embedded star-forming cores 
with masses M core , the radius r s , the mean gas number 
density n#, and the r.m.s velocity dispersion v rms of the 
cores, as well as the final stellar mass m*/=e core M core , 
can be estimated based on the following formulae (details 
see McKee & Tan 2002): 



n H 



, =O.O74(m* / /3OM ) i ^E- i /^ pc, 
1.0 x lO^m^/SOMo)- 1 / 2 ^ 3 / 2 en 



(5) 
(6) 



v rrns = 1.65(m* f /30M Q ) 1/4 Z 1/4 kms" 1 , (7) 

where E is the mean clump surface mass density in units 
of gem -2 , and e core (assumed to be 0.5 here) is the frac- 
tion of the core mass which is eventually accreted onto 
the central star. From Figure [TJ the majority of the 
densest clumps are embedded in the region above the 
12cr contour level (i.e. 288 mJybeam -1 ), which corre- 
sponds to the E of 0.31-0.76 gem -2 . Based on the val- 
ues in Table [TJ we adopt a fiducial value n#~10 6 cm -3 
for an example. We estimate the core radius r s ~0.023- 
0.056 pc (0.58"-lA"), the r.m.s. velocity dispersion 0.51- 
1.3 kms -1 , and the final stellar mass m*f~0. 89-13 M . 
The former two quantities can be examined with <0.5" 
resolution Atacama Large Millimeter/submillimeter Ar- 
ray (ALMA) observations in the near future. The cor- 
responding core mass M core =2xm*/ in our estimates is 
several times smaller than the masses of the parent gas 
clumps listed in Table Q] so may be reasonable. 

As discussed in the previous section, the dense clumps 
embedded in some regions (e.g. A, B, D, E, F, G) are as- 
sociated with the 22 GHz water maser and the 36.2 GHz 
and 44.1 GHz Class I CH3OH masers, which are often 
seen in the early phase of high-mass star-formation. In 
addition, the earlier VLA observations of the 2 cm and 
the 6 cm continuum emission have found ultracompact 
Hll regions embedded with several OB stars in the east 
of the 50 kms -1 cloud (e.g. Ho et al. 1985). The ob- 
served high HCO + 4-3 velocity dispersion may be inter- 
preted by (proto) stellar feedback. However, Sjouwerman 
et al. (2010) also suggested that the 36.2 GHz Class I 
CH3OH masers can be excited at the post shock regions 
created by cloud-cloud collisions. High angular resolu- 
tion molecular line observations and studies of the maser 
proper motion may distinguish these two cases, although 
they are not mutually exclusive. In fact, it has been 
argued that the predominant mode of star formation is 
via external compression of molecular clouds by cloud 
collisions and supernova (Morris 1993). The high-mass 
star formation can also be induced by AGN activities 
(Silk et al. 2012). In these cases, the physical proper- 
ties of the molecular gas reservoir feeding the high-mass 
star formation may be very different from the typical OB 
star- forming cores. 
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Fig. 10. The velocity integrated SiO/C 34 S 1-0 ratio and the 
0.86 mm intensity at pixels where the velocity integrated C 34 S 1-0 
intensities are higher than 0.7 Kkms -1 (gray dots). The values at 
the locations of the 1720 MHz OH masers (green), the 1612 MHz 
OH masers (light blue), the 22 GHz water masers (black), the 36.2 
GHz CH3OH masers (magenta), and the 44.1 GHz CH3OH masers 
(blue) are presented by colored symbols. For maser, the lower limit 
of the velocity integrated SiO/C 34 S 1-0 ratio are given if the C 34 S 
1-0 emission is not robustly detected. The data points located 
outside of the field of view of either the 0.86 mm image or the SiO 
image are not presented. 



4.4. The Non- Uniform CND 

Previous interferometric spectral line observations 
have shown abundant localized structures in the cen- 
tral 2-4 pc CND (see Christopher et al. 2005; Montero- 
Castaho et al. 2009; Martin et al. 2012, and references 
therein). The distributions of the spectral line emis- 
sion, however, are sensitive to the local temperature, 
volume density, chemical abundances, and other radia- 
tion transfer effects such as the foreground absorption or 
self-absorption. The optically thin 0.86 mm dust ther- 
mal continuum emission is a more robust tracer of the 
gas mass (e.g. Hildebrand 1983). Without being sub- 
jected to missing flux, our JCMT+SMA 0.86 mm dust 
continuum image (Figure [T]) successfully reproduces the 
clumpy CND structures seen in the spectral line obser- 
vations (e.g. Figure HJ . It appears that some previously 
known clumps in the CND are not merely due to ex- 
citation or radiation transfer effects. For example, the 
0.86 mm emission clump A and Bn (Table [TJ Figure [3]) 
coincide with the HCN 4-3 emission clump CC and BB 
reported by Montero-Castaho et al.(2009). As can be 
expected, we found that around the CND, the velocity 
integrated HCO + 4-3 intensity is correlated with the in- 
tensity of the 0.86 mm continuum emission (Figure fTTj) . 
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Fig. 11. — The velocity integrated HCO+ 4-3 intensity and the 
0.86 mm intensity of pixels within the 60" radius around the Sgr 
A*. We color code the pixels in the four selected zones in Figure 
[T2l We note that our HCO + 4-3 image is subjected to missing flux, 
however, should not significantly bias the analysis of the emission 
from the localized clumps. 



However, in Figure [TTJ we see that the correlation be- 
tween the HCO + 4-3 emission and the 0.86 mm emission 
is dominated by at least two distinct populations which 
can be discerned with our current angular resolution and 
sensitivity. We selected the high HCO + 4-3 intensity 
parts of the two dominant populations as Zone 1 and 
Zone 2 in Figure [TTJ and additionally one high HCO + 4-3 
intensity and high 0.86 mm intensity population as Zone 
3, and one high 0.86 mm intensity population as Zone 4. 
We argue that these high intensity pixels are more likely 
to be associated with compact bright clumps, for which 
the HCO + 4-3 intensity is less biased by missing flux. 

The spatial distributions of the pixels in these four 
zones are not random. For example, the pixels in Zone 
1 seems to trace the two arc-shaped clumps west of the 
CND, the eastern edge of the CND, and the clumps lo- 
cated northwest of the Sgr A*. The pixels in Zone 2 are 
associated with the Northeast Lobe, the Southern Ex- 
tension, and northern part of the Southwest Lobe. The 
southern part of the Southwest Lobe is associated with 
the pixels in Zone 4. Comparing Figure [12] with the ratio 
of the HCN 4-3 to HCN 1-0 integrated intensity reported 
by Montero-Castaho et al. (2009), we deduce that the 
very dense clump in the southern part of the Southwest 
Lobe has a low gas excitation temperature. This provide 
hints for the the cooler exterior gas clumps raining down 
on the warmer CND. For a dynamically evolving CND 
which is connected with several exterior gas streamers 
(see Liu et al. 2012), it is not surprising that the excita- 
tion conditions, the chemically abundances, or the dust 
properties are not yet homogenized. 

5. SUMMARY 

We present a wide-field SMA mosaic toward the Galac- 
tic center. We also observed the CS/C 34 S/ 13 CS 1-0 and 
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the SiO 1-0 lines using the GBT. The optically thin 0.86 
mm dust thermal continuum image with ~5" angular res- 
olution confirms the 2-10 pc scale gas streamers and the 
detailed structures of the CND, which were seen in pre- 
vious molecular line observations. We marginally resolve 
more than 22 massive (10 1 -10 3 M ) gas clumps, which 
are embedded in very turbulent molecular gas clouds 
(Sv ~ 10-30 kms -1 ). Examination of the brightness ratio 
of the HCO + 4-3 line and the 0.86 mm continuum emis- 
sion shows that the non-homogenized central 2-4 pc CND 
has dense clumps with distinct excitation conditions or 
chemical abundances. While the distributions of the 1720 
MHz and 1612 MHz OH maser clusters do not show ob- 
vious correlations with the dust emission, the 22 GHz 
water masers and especially the 36.2 and 44.1 GHz Class 
I CH3OH masers are seen preferentially near the dense 
clumps. Even the most significant dense clumps in our 
selected samples are marginally unstable against the tidal 
force, and we do not find any self-gravitationally bound 
gas clump associated with the Class I CH3OH masers or 
the 22 GHz H2O masers. How the OB stars form in such 
environment remains puzzling. The mechanisms to form 
the dense clumps and the high-mass stars in the Galac- 
tic center might be very different from the mechanisms 
in typical giant molecular clouds. Our simple estimates 
suggest that if the detected dense clumps are confined 
by the high external pressure, the presumably existing 
embedded virialized gas cores can form high- mass stars. 
Deep JVLA observations to search for the signature of 
the stellar photoionization, or ALMA observations of the 
high excitation hot core tracers to look for the gravita- 
tionally accelerated rotation, may diagnose the OB star 
formation in the dense clumps. 

Our GBT data suggest mildly enhanced SiO/C 34 S line 
ratio towards the Southern Arc and the 20 kms -1 cloud, 
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relative to the SiO/C 34 S line ratios in the other gas 
streamers in the observed field. Our observations do 
not yet resolve the recognizable correlation between the 
SiO/C 34 S line ratio and the intensity of the submillime- 
ter dust emission. Higher angular resolution and more 
sensitive observations are required to examine whether 
the formation of the dense gas structures are predomi- 
nantly induced by shocks. 

The GBT and SMA data are from projects 
GBT11B050, SMA2011AS085 and SMA2011BS040, 
which are parts of the integrated state-of-art imaging 
project KISS: Kinematic Processes of the Extremely 
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APPENDIX 
A. SHORT SPACING INFORMATION 

Short spacing information, below ~7.0 &A, was obtained from single dish telescope observations. We retrieved the 
archival public processed JCMT0 SCUBA (ProjectID: M01BU26, observed on 2001 August 05) 0.86 mm continuum 
image (Figure [13j rms~50 mJybeam -1 ). We used the MIRIAD task DEMOS to generate SMA primary beam weighted 
models for each SMA pointing. For each SMA field, we then resampled the primary beam weighted models in the 
uv domain with 185 Gaussian randomly distributed visibilities using the MIRIAD tasks UVRANDOM and UVMODEL. We 
manually assigned a system temperature of 350 K to the visibility model to adjust the weighting. 

The typical pointing accuracy for JCMT SCUBA is about 3" and the tracking accuracy ~1.5" (Di Francesco et al. 
2008). From their JCMT SCUBA legacy survey, Di Francesco et al. (2008) also mentioned that larger pointing offsets 
(e.g. ~6") occurred occasionally. When combining with the interferometric observations, pointing offsets during the 
single dish observations can cause defects especially near bright compact sources. We used a limited uv range of 0-3.8 
kA for the single dish uv model to suppress the potential effect of the single dish pointing offsets, as well as the effect 
of the single dish primary beam and deconvolution errors. Although the effect of the single dish primary beam is not 
fully eliminated, it is smaller than the ~20% absolute flux error in typical SMA observations (Appendix [B]). Limiting 
the uv range implies that the single dish primary beam causes defects primarily for structures >33", and therefore 
should not confuse the identification of local gas clumps in the streamers. 




40° 
RA (J2000) 

Fig. 13. — The JCMT SCUBA 0.86 mm continuum image (0 ma j x <9 min =14 // .3xl4 // .3). Contour spacings are 5a starting at 5a (a=50 

mJybeam - 1 ). For our convenience in discussion, we mark the peaks above the 45a significance level besides the Northeast Lobe and 
Southwest Lobe as JCMT-1,2,3. The symbols are described in Figure [T] 

We jointly deconvolved and imaged the concatenated single dish uv model and the SMA data using the MIRIAD 
software package (using tasks: INVERT, M0SSDI, REST0R). The gap (3.8-7 kA) in the uv sampling may yield uncertainties 
in reconstructing the 14 // -33 // angular scale structures. However, structures with these angular scales can be directly 
inspected from the single dish image. We generated the high angular resolution continuum image with the parameters 
robust=0 fwhm=4,4 m INVERT, which yield a synthesized beamwidth ^ 



ma] 



5".lx4".2, and BPA. -28°. To 
compare with the earlier observations, we also generated one lower angular resolution (# m aj x # m in=7 // -4x6 // .2, BPA. 
~44°) continuum image with the weighting parameters robust=0 fwhm=8,8 to enhance the sensitivity to the extended 
emission. 



3 The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Facil- 



ities Council of the United Kingdom, the Netherlands Organisa- 
tion for Scientific Research, and the National Research Council of 
Canada. 



15 



Jy beam 1 
-0.4-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



Jy beam 
-0.4-0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 




17 h 45 m 52 s 48 s 



40 
RA (J2000) 



44° 40 s 36 s 

RA (J2000) 



Fig. 14. — Left: the residual image constructed by subtracting the JCMT SCUBA image from the smoothed high angular resolution 
SMA+JCMT 0.86 mm continuum image (see Appendix [B]) . The JCMT beam is shown in the lower left. Contours are 3crx[-2, -1, 1, 2, 3, 
4, 5] (cr=89 mJybeam -1 ). We note that the averaged 3a flux density 267 mJybeam -1 in this image corresponds to the averaged flux of 
~28 mJy in each 5 // .lx4 // .2 SMA synthesized beam. The symbols are described in Figure[T] Right: similar to the left panel. However, the 
original JCMT SCUBA image was shifted toward the south by A" before being subtracted from the SMA+JCMT image. 



B. CONSISTENCY CHECK AND NOISE STATISTICS 



We smoothed the high angular resolution SMA+JCMT image to the angular resolution of JCMT SCUBA, and 
then subtracted the JCMT SCUBA image from the smoothed SMA+JCMT image. The residual image after the 
subtraction is shown in Figure [HJ The lower angular resolution SMA+JCMT image, and the previously published 
SMA HCN 4-3 image (Liu et al. 2012) are presented in Figure [15] On average the contour intervals in Figure 
[141 (i.e. 267 mJybeam -1 ) correspond to the la (24 mJybeam -1 ) detection level of the higher angular resolution 



(0, 



maj 



#min=5 // .lx4 // .2) SMA+JCMT 0.86 mm continuum image (see below). 



Our overall impression by comparing Figure HH and [15] is that the limited uv sampling rate of our SMA observations 
cause imaging defects near the bright compact sources. This is most obviously seen around the Sgr A*, the Southwest 
Lobe, and the bright clump in the Southern Arc (coincide with JCMT-1). However, we found that manually shifting 
the JCMT image towards the south by 4" can suppress the residual (e.g. Figure HH]) around the Sgr A*. In our zero 
spacing model (see Section [Aj), we did not manually correct for this 4" offset for the sake of objectiveness. We do 
not think this potential JCMT 4" pointing offset can cause very obvious imaging defects since only the >33" scale 
structures are taken for modeling the zero spacing. However, it can lead to the mismatch between the combined 
SMA+JCMT image and the JCMT image, which will result in some residual in Figure fl4l 

Beyond the aforementioned regions, the extended structures appear to be reasonably reconstructed in the 
SMA+JCMT image, and the residual image has an rms noise level of ~89 mJybeam -1 . The Northern/Southern 
Ridges, Northeast/Southwest Lobes, and the W-2,3,4 Streamers are detected in the lower angular resolution 
SMA+JCMT image (Figure [15] left). The extended residual features at the edge of Figure [14] are caused by the 
sidelobe features in the JCMT SCUBA image (e.g. the negative brightness sidelobe at northwest of Figure [T3l causes 
the positive excess in Figure [T4|) . These single dish sidelobes are extended (>15") and should not confuse the imaging 
of localized gas structures. However, locally they can bias the absolute brightness in the higher angular resolution 
SMA+JCMT image (mostly only 2a). The >33" scale defects caused by the single dish primary beam effect discussed 
in the previous section appear lower than the la detection level of the SMA+JCMT image. 

The theoretical rms noise level of the SMA+JCMT image is ~14 mJybeam -1 . Because of the absence of emission- 
free areas in our map, it is difficult to directly measure the rms noise level we actually achieved. Since only the >33" 
structures in the JCMT image, which have high signal to noise ratios, are combined with the SMA data. We assumed 
the noise in the SMA+JCMT image only has a weak dependence on the noise in the JCMT image. We therefore can 
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Fig. 15. — The combined SMA+JCMT 0.86 mm continuum image (color). Left: Color image shows the tappered SMA+JCMT 0.86 mm 
continuum image (# ma i X # m j n =7 // .4x6 // .2). The synthesized beam of the SMA+JCMT image is shown in the lower left. Contours show 

the velocity integrated SMA HCN 4-3 image taken from Liu et al. (2012) (# ma j x # m i n =5 // .9x4 // .4). Contours are 75 Jy beam _1 kms _1 x [1, 
2, 4, 8, 16]. By comparing with the GBT CS 1-0 spectra, Liu et al. (2012) suggested that this HCN 4-3 image is subjected to the significant 
missing flux issue around the 50 km -1 cloud, the Molecular Ridge, the 20 kms -1 , and the Southern Arc. Right: Color image shows 
the high angular resolution SMA+JCMT 0.86 mm continuum image (# ma i x # m j n =5 // .lx4 // .2). Contours show the free-free continuum 



emission model (see Appendix [C| , which is smoothed to the same angular resolution with the SMA+JCMT image. 
mJybeam- 1 x[l, 2, 4, 8, 16, 32, 64]. 

measure the rms noise level of the high angular resolution SMA+JCMT image by 

\f (rms of the residual imagej — (rms of the JCMT imagej 
y (JCMT beam area)/(SMA+JCMT synthesized beam area) 



Contours are 10 



(Bl) 



This yields an rms noise level of 24 mJybeam -1 (0.011 K) in the high angular resolution SMA+JCMT image, which 
is ~1.7 times higher than the theoretical noise level. Generally, the SMA+JCMT image is still dynamic-range limited, 
especially near the bright compact objects. This effect hampers the systematical search and statistical study of the 
clumpy structures. 

C. FREE-FREE MODEL SUBTRACTION 

The emission from the ionized mini-spiral arms can be recognized from the high angular resolution SMA+JCMT 
image (Figure [T5|). We retrieved the archival VLA Q band observations taken on 2003 February 14, March 17, April 
16, and 2004 March 14, April 21, May 16 to generate a model of 0.86 mm free- free emission. These VLA observations 
were taken with dual 50 MHz IFs centered at 43.314 GHz and 43.364 GHz, with full polarization. The sampling range 
of these VLA data, is 3.7-490 kA. However, we assume that only the localized bright components will significantly 
contribute to the 0.86 mm free- free continuum emission. The basic calibrations and self-calibrations were performed 
using the Astronomical Image Processing System (AIPS) software package of NRAO. We scaled the VLA image such 
that the central point source Sgr A* has the flux of 1 Jy (Bower & Backer 1998; Yusef-Zadeh et al. 2011), and then 
smoothed the VLA image to the angular resolution of the SMA+JCMT image. The smoothed and rescaled VLA 
image, as the zeroth order model of the 0.86 mm free-free continuum emission, is presented in the right panel of Figure 
[T5l The brightness of the 3 mm and the 1.3 mm emission is comparable, which can be checked in Kunneriath et al. 
(2012). We then subtracted the free- free emission model from the SMA+JCMT image (Figure [1]). As can be seen from 
Figure [15] and [TJ this free-free model subtraction can only manifestly change the geometry in the Northern and the 
Eastern ionized mini-spiral arm regions and around Sgr A* (see Zhao et al. 2009 and references therein for the ionized 
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mini-spiral arm). We note that the scaling of the VLA image in this process cannot take into the consideration of the 
spatial variation of the spectral index (Kunneriath et al. 2012). We have checked that if we scaled the VLA image 
to be >20% brighter, this free-free model subtraction will induce noticeable over-subtracted features in the Northern 
and the Eastern ionized mini-spiral arms. Except for the location of Sgr A*, in the central V area, the free- free model 
subtracted SMA+JCMT image provides a lower limit for the 0.86 mm dust thermal emission. 



